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ZEROS OF HARMONIC POLYNOMIALS, CRITICAL LEMNISCATES AND 

CAUSTICS 

DMITRY KHAVINSON, SEUNG-YEOP LEE, ANDRES SAEZ 


Abstract. In this paper we sharpen significantly several known estimates on the maximal 
number of zeros of complex harmonic polynomials. We also study the relation between the 
curvature of critical lemniscates and its impact on geometry of caustics and the number of 
zeros of harmonic polynomials. 
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1. Introduction 

We concern ourselves in this paper with complex harmonic polynomials , i.e., polynomials 
which admit a decomposition 

h{z) = p{z ) +q(z), 

where p = p n and q = q m are analytic polynomials of degrees n and m respectively. An interesting 
open question is, for given n and m, to find the maximal number of solutions to the equation 
h{z) = 0, i.e., extending the Fundamental Theorem of Algebra to harmonic polynomials, see [8] 
and references therein. Throughout the paper we assume n > to, for the case n = m could give 
an infinite solution set. 

Wilmshurst [17, 18], in his doctoral thesis, proved the following: 

Theorem (Wilmshurst). The equation h(z) = 0 has at most n 2 solutions. 

The proof of this result readily follows from Bezout’s theorem [3]. 

Seeking to improve on this bound, Wilmshurst conjectured that the equation h(z) = 0 has 
at most 3n — 2 + m(jn — 1) solutions. Khavinson and Swiqtek [ ] confirmed Wilmshurst’s 
conjecture when m = 1 using complex dynamics, and the bound was shown to be sharp by 
Geyer [5]. However, as was shown by Lee, Lerario, and Lundberg [)], the conjecture is not true 
in general, for example, when m = n — 3. Also see [6] for many more counterexamples. 

Our first theorem bounds the number of roots off the coordinate axes for the harmonic poly¬ 
nomials with real coefficients. 


l 
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Theorem 1. For a harmonic polynomial h(z) = p n (z)+q m (z) with real coefficients, the equation 
h(z) = 0 has at most n 2 — n solutions that satisty (Re z) (Iinz) ^ 0. 

Our next two theorems provide lower bounds on the maximal number of roots. 

Theorem 2. For all n > m, there exists a harmonic polynomial h{z) = p n (z ) + q m (z) with at 
least 3n — 2 roots. 

Theorem 3. For all n > m, there exists a harmonic polynomial h{z) = p n (z) + q m (z) with at 
least to 2 + m + n roots. 

The above three theorems will be proved in Section 2. 

Remark. 

(i) The reason why in Theorem 1 we only consider roots off the coordinate axes is the 
following. Consider pfz) = z n + (z — l) n and q(z) = z n — (z — l) n . Then h{z) = p(z)+q(z) 
has n 2 number of roots including the root at 0 with the multiplicity n. In fact, this is 
the polynomial that Wilmshurst used (with a slight perturbation to split the multiple 
root at the origin) to show that the maximal bound n 2 is sharp. This example shows 
that Theorem 1 is sharp. Theorem 1 also yields that the harmonic polynomial with real 
coefficients and with the maximal ( n 2 ) number of roots should have at least n roots on 
the coordinate axes as in the above example. 

(ii) Theorem 2 and Theorem 3 complement each other. Theorem 2 is stronger than Theorem 
3 when m 2 +m+2 < 2 n and Theorem 3 is stronger than Theorem 2 when m 2 +m+2 > 2 n. 
Also, Theorem 2 is not at all trivial since the argument principle for harmonic functions 
[4, 14, 16] only yields that h has at least n zeroes (see Fact 1 in Section 2.2). 

(iii) Note that, compared to Wilmshurst’s conjecture, Theorem 3 undercounts the number of 
roots by 2 (n — m — 1). See Section 2.4 for the in-depth discussion. 

Theorem 3 allows the following important corollary. 

Corollary 1. Let Z ntTn denote the maximal possible number of zeros of h = p n + qfd- Then, for 
any fixed integer a > 1, we have 

Z nn — a 

hmsup-^— = 1. 

n—> oo Tl 

More generally, if m = an + o(n) with 0 < a < 1, we have 

v Z n m 2 

lim sup —— > a . 

n—f oo n 

Proof. From Theorem 3 we have Z n ^ m > n 2 + 2(1 — a)n + a(a — 1) for the first case, and 
Z n ,m > a 2 n 2 + o(n 2 ) for the second case. ■ 

This corollary yields that the maximal number of roots is asymptotically given by Wilmshurst’s 
theorem. This answers a question posed by the first author more than a decade ago. Also this 
corollary complements the estimates on the expected number of zeros of Gaussian random har¬ 
monic polynomials obtained by Li and Wei in [11] and, more recently, by Lerario and Lundberg 
in [10]. For example, for m = an, a < 1, the expected number of zeros is ~ n ([LI]) for the 
Gaussian harmonic polynomials and ~ c a n 3 / 2 [10] for the truncated Gaussian harmonic polyno¬ 
mials. Yet, Corollary 1 yields that among all harmonic polynomials the maximal number ~ an 2 
of zeros occurs with positive probability, thus expanding further the results in [2] for m = 1. 

To state our last theorem, we have to introduce the set 

^ = {z : |/(z)| < 1}, 


(1) 
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where 


/(-) = 


KO ) 


Recall that the mapping z i—> h(z) is sense-reversing precisely on fi, i.e. the Jacobian of the map 
h is negative on ft. The boundary dfl is the lemniscate {z : |/(z)| = 1}. 

Each connected component of f 1 must contain at least one critical point of p n . Indeed, if there 
is a connected component without a critical point, by applying the maximum modulus principle 
to f(z) and 1 /f{z) with \f(z)\ = 1 on the boundary of that component, we have that / is a 
unimodular constant, a contradiction. 

This implies that there are at most degp' n = n — 1 connected components of ft. 

For m = 1 (when Wilmshurst’s conjecture was proven to hold [ f]), Wilmshurst guessed ([17], 
p.73) that the following might be true: “In each component of ft where h(z) = p n {z) +z is sense 
reversing , the behavior of h will be essentially determined by the z term so there will only be one 
zero of h” . From that the maximal number of roots (i.e., 3n — 2 for m = 1) may be obtained by 
the argument principle when each connected component of ft contains a zero. 

The statement above is true when the component of f 1 is convex according to the following 
result (see [17], p.75). 


Theorem (Sheil-Small). If g(z) is an analytic function in a convex domain D and |</( 2 )| < 1 
in D, then z + g(z) is injective on D. 


Note that the theorem relates the geometry of critical lemniscates with the number of zeros, 
because D can have at most one zero of the function ~z + g(z) if the latter is injective. In a non- 
convex component of ft, it is possible to have more than one zero of h(z). In [17] an example 
is given where a non-convex component of ft contains two critical points of p n and two zeros 
of h{z) = p n {z ) + z. This is not surprising because each component of f l can have a zero of h 
and, therefore, one can have two zeros in a component by merging two components into one. 
The resulting component then has two critical points of p n . It was however not clear whether a 
connected component containing a single critical point of p n could possibly have more than one 
zero of h. 

Here we show that it is indeed possible and, also, present a necessary and sufficient condition 
for having more than one zero of h in a connected component of ft that contains a single zero of 
/ (i.e., a single critical point of p n ). The theorem holds for general m and n. 

Theorem 4. Let n > m. Let D be a connected component of It (defined by (l)j containing 
exactly one zero of f. On a smooth part of the curve q m (dD) (the image of dD under q m ) let 
k be the curvature of q m (dD) with respect to the counterclockwise arclength parametrization of 
dD. Then, the following are equivalent: 

i) Let f{z) = p' n (z) / q' m (z). There exists z G dD such that 

k(z) ^ 1 

ii) There exists 9 £ R. and A € C such that 

Pn{z)~ q m {z) 

has at least two zeros in D, where p n {z) = e te p n (z) + A. 


Remark. 

(a) Note that k is the curvature of dD when q m (z ) = z. In this case, the theorem tells us 
exactly how “non-convex” the domain D needs to be in order to have multiple zeros of 
h, improving upon Sheil-Small’s theorem. 
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(b) In statement (ii) of Theorem 4, we note that |pv>(z)| = \p' n (z)\. The corresponding 
lemniscate, {z : \j? n {z) / q' m {z)\ = 1}, is therefore, the same for all 9 and A. 

For (n,m) = (4,1), we provide an example where a component of Q with one critical point 
of p n contains two zeros of h, see Figure 1. Note that the roots appear where is (slightly) 
concave. The example is produced based on the discussion following our final Theorem 5 in 
Section 4 regarding the shapes of critical lemniscates. 



Figure 1. Roots (dots) for h(z) = -(0.934124 + 0.356949*)+* + (0.0581623 + 
0.156514f) 2 4 + (0.354765 - 0.131835*)z 3 - (0.116325 + 0.313028*)z 2 - (1.06429 - 
0.395504*)2 + (0.247627 + 0.020994*). The stars are for the zeros of / and the 
shaded region is ft. The right picture is a zooming image of the top component 
of Q . 


Acknowledgement. This work resulted from an REU group discussions that also included 
Prof. Catherine Beneteau and Brian Jackson. The second author was supported by Simons 
Collaboration Grants for Mathematicians. The first and third authors were supported by the 
USF Proposal Enhancement Grant No. 18326 (2015), Pis: D. Khavinson and R. Teodorescu. 
We are greatly indebted to the referee for pointing out an error in the initial version of Theorem 
1 that was due to a misinterpretation of D. Bernstein’s theorem. 

2. Proofs of Theorems 1, 2 and 3 
2.1. Proof of Theorem 1. Given a bivariate real polynomial P defined by 

n m 

P{x, y) = a t 7 X V) a ij e 

i—0 j—0 

the Newton polygon Afp of P is the convex hull of Np C M 2 where Np = {( i,j ) : ay ^ 0}. 

Given two polynomials P and Q with Newton polygons A fp and ATq, let Mp,q be the 
Minkowski sum of Afp and AAq , defined by 

Mp,Q = {(*1 + * 2 ) ji + J 2 ) I (h; Ji) € Afp, (* 2 ? J 2 ) £ A/q}. 
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Figure 2. The Newton polygons A (left) and B (right) for n = 6 (top) and 
n = 5 (bottom), respectively. 


Let [X] denote the area of a set X C M 2 . We then define the mixed area of P and Q as 
[Mp,q\ - [Up] - [ATq]. 

Theorem (D. Bernstein; cf. [15], or the original articles [1, 13]). The number of solutions to 
the system of polynomial equations p(x, y) = q{x , y) = 0 satisfying xy ^ 0 does not exceed the 
mixed area of p and q. 

For any analytic polynomials p n {z) and q m (z) of degree n,m, we have 

h{x + iy) =p n (x + iy) + q m {x + iy) = A(x, y) + iB(x, y), 

where A and B are polynomials with real coefficients. Let A be the Newton polygon of A(x, y) 
and B be the Newton polygon of B[x , y). As we will see below, the polynomials p n and q n having 
real coefficients leads to a certain structure of A and B. 

Lemma 1. Given a generic h(z) with only real coefficients let A and B be defined as above. If n 
is even, A is the isosceles triangle with vertex set {(0, 0), (0, n), (n, 0)} and B is the trapezoid with 
vertex set {(0,1), (0, n — 1), (1, n — 1), (n — 1,1)}. If n is odd, A is the trapezoid with vertex set 
{(0,0), (0, n —1), (1, n— 1), (n, 0)} and B is the isosceles triangle with vertex set {(0,1), (0, n), (n— 
1,1)} (see Figure 2). 


Proof. The lemma follows directly from the fact that 

A(x,y) = ^2 aj k x 3 y k , B(x,y) = ^ b jk x J y k , a jk ,bj k £R. 

j-\-k<n, j-\-k<.n, 

k is even k is odd 

The condition, j + k < n, on the summation indices follows from deg h = degp„ = n. The 
condition on k to be even or odd comes from p n and q m having real coefficients. (Note that in 
the expansion of p n (x + iy) or, of q m {x + iy), the “i” comes only together with y - for a term 
of the form x 3 y k , the power of “i” is given by k, the power of y.) 
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We need to show that all the coefficients that correspond to the extreme points of the convex 
hulls are nonvanishing. First of all, those that satisfy j + k = n: 

ao, n , a n> o, &i, n -i and b n -i,i for even n, 
ai,n-i, a n ,o, b 0 , n and b n _ 1A for odd n, 


are all nonvanishing because deg p n = n. 
Using h(z) = p n (z ) + 9m(~) 5 we note that 


h(x), 


d_ 

dy 


h(x + iy) 


v —o 


d n ~ x 

g^i h (x + iy) 


are all nontrivial polynomials in x, due to the presence of the terms, x n , x n ~ 1 y and xy n ~ x in 
h(x + iy) respectively. So there exists ij € R such that, when x = Xq, none of the above three 
polynomials vanishes. Defining 


h(z) = h(z + x o), xq € R, 


h. is a real polynomial (of the same holomorphic and antiholomorphic degrees) and has the same 
number of zeros as h. Moreover, the coefficients at the extreme points are all nonvanishing 
because 


a 0 ,o = h( 0) = h(x 0 ), 


d ~ 

bo,i = w~h(iy) 
dy 


d 


Vo,n-l — 


1 


(n- 1)! dy 


= w-/i( X 0 + iy) 

v =o Q y 

d n ~ x ~ 
zi h(iy) 


y =o 

1 


y=0 


(n— 1)! dy 


d n -i 

— h(x 0 + iy) 


y=0 


where, in the last line, the symbol stands for a, when n is odd, and for b , when n is even. 


Now we prove Theorem 1. We divide the proof into two cases. 

Case 1. Let n be even. From Lemma 1, A is an isosceles triangle with vertex set {(0, 0), (0, n), (n, 0)} 
and B is the trapezoid with vertex set {(0,1), (0, n— 1), (1, n— 1), (n— 1,1)}. Now, the Minkowski 
sum M of A and B is a trapezoid with vertex set {(0,1), (0, 2 n— 1), (1, 2 n— 1), (2n— 1,1)}. The 
mixed area is, then, 

[M] - [A] [B] = (2n 2 - 2n) - y - 

Case 2. Let n be odd. From the lemma, A is the trapezoid with vertex set {(0,0), (0, — 

1), (1, n — 1), (n, 0)} and B is the isosceles triangle with vertex set {(0,1), (0, n), (n — 1,1)}. The 
Minkowski sum A4 of A and B is then the trapezoid with vertex set {(0,1), (0,2 n — 1), (1, 2 n — 

1), (2 n — 1,1)}. The mixed area is 

[M\ - [A] [B] = (2n 2 - 2 n) - - (? ^ )2 


2.2. Proof of Theorem 2. For a harmonic function h(z) 
sense-preserving at a point 2 if the Jacobian of h, 


det 


c^Re h{x + iy) 
clr Im h(x + iy) 


d y Re h{x + iy) 
dylm h(x + iy) 


x-\-iy=z 


= p{z) + q(z), we say that h is 


= \p\z)\ 2 -\q'{z)\\ 


is positive, and sense-reversing at z if h(z) is sense-preserving at z. Otherwise, we say h is 
singular at z. We also say that h is regular if all the zeros of h are either sense-reversing, or 
sense-preserving. 
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For an oriented, closed curve T such that a continuous function F does not vanish on F, we 
denote by Ap arg F(z) the increment in the argument of F along T. The following is well known. 

Theorem (The argument principle for harmonic functions [1, 14]). Let FI be a harmonic function 
in a Jordan domain D with boundary T. Suppose FI is continuous in D and H ^ 0 onF. Suppose 
FI has no singular zeros in D, and let N = N + — N_, where N + and N— are the number of sense¬ 
preserving zeros and sense-reversing zeros of H in D respectively. Then, Ap arg H{z) = 2 ttN. 

The next fact [ ] follows then by applying the argument principle to a circle of a sufficiently 
large radius where \p n \ \q m \- 


Fact 1. Let h = p n + q m be regular. Let N + be the number of sense-preserving zeros of h and 
N _ be the number of sense-reversing zeros. Then, 

n = N + — N_. 

An elegant proof (due to Donald Sarason) of the following lemma can be found in [7]. 

Lemma 2. If p{z) is a polynomial of degree greater than 1, then the set of complex numbers c 
for which p(z) + q(z) — c is regular is open and dense in C. 

For m = 1, Bleher et al. [ ] proved that in the space C n+1 of harmonic polynomials p n (z) + z 
of degree n > 2 the set of “simple polynomials”, i.e., regular polynomials with k zeros, is a 
non-empty open subset of C n+1 if and only if k = n, n + 2,..., 3n — 4,3n — 2, and that the 
non-simple polynomials are contained in a real algebraic subset of C n+1 . 

The following fact is noted in [12], Theorem 3.2. 

Lemma 3. If the function h(z) = p n {z) + z has 3n — 2 zeros, then h(z) is regular. 


Proof. Assume h is not regular and has exactly 3n — 2 roots. The number of roots that are not 
sense-preserving is at most n — 1 because each of those roots attracts a critical point of p n under 
the iteration of z ha p n (z), cf. Proposition 1 in [ ]. Therefore, the number of sense-preserving 
roots must be at least 2n — 1. For all non-singular roots zfs, there exists e > 0 such that the 
disks, B t (zj) = {z : \z — zj\ < e}, do not intersect each other, do not intersect the singular set 
{z ■ \p'{z) | = 1}, and A dB,( Zj ) arg h = ±2t r. Defining 


6 = 


min 

z£dB e (zj) 


\Kz)l 


for any c £ C with |c| < S, we have A^s E ( Zi ) arg h = A dB c ( Zj ) arg {h — c) and, using the argument 
principle for harmonic functions, the equation h(z ) — c has exactly one zero in each B e (zj). 

Suppose zq is a singular zero of h. Since h(zo) = 0 and h is continuous near zq, there is an 
q > 0 such that B v (zq) does not intersect any B e (zj) and 

\h{z)\ < S for all z £ B v (zq). 

Further, the set B ri (z 0 ) intersects sense-preserving region, since otherwise, log|p'(z)| would be 
constant over B v (zq) by the maximum modulus principle. Let f be a sense-preserving point in 
B v (zo). We can set c = h(f) since \h(()\ < S, and h(z) — c = h(z) — h(f) must have zeros in 
each B e (zj) and at f £ B v (zq). Consequently, h(z) — h{Q has at least 2 n sense-preserving zeros. 
By Lemma 2, we can choose ( such that h(z) — h(f) is a regular polynomial, which contradicts 
the result of Khavinson and Swiqtek [7] that the regular polynomial can have at most 2 n — 1 
sense-preserving roots. ■ 


We now complete the proof of Theorem 2. 
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Proof. Let p n (z) be an analytic polynomial such that the equation p n (z) + z = 0 has 3n — 2 
solutions. By Lemma 3 the polynomial h(z ) = p n (z) + z is regular. Let Zq be a zero of h. 
One can find a circle, T, centered at zq with radius e such that h does not vanish on T and 
Ar arg h = ±27r. As in a standard proof of Rouche’s theorem, taking S such that 

max zer (|z| m + 1) 

we have the perturbed mapping, z i—> h(z) + <5z m , that preserves the winding number of h(T), 
that is, the perturbed mapping also vanishes in the interior of T. Applying the same argument 
to all 3 n — 2 zeros of h , we can choose <5 such that h{z) + 6z m has (at least) 3n — 2 zeros. Setting 
now q m (z) = Sz m + z completes the proof. ■ 

Remark. This proof suggests that if the equation p n (z) + g m _i(z) = 0 has at most k solutions, 
then there exists a harmonic polynomial p n (z) + q m (z) with k zeros. However, a proof of this, 
following the above argument, would require that p n {z) + q m -i{z) be regular. 

2.3. Proof of Theorem 3. Let us sketch the procedure (similar to [9, 6]) that allows creating 
examples of harmonic polynomials with a large number of roots (cf. Figure 3). 

Fix n and m < n. Let 

S(z) = (z — a) ra-1 (z + (n — l)a), 

T{z) = (z- fer+V-™- 1 + t n . m . 2 z n - m - 2 + ■■■+ to). 

The n — m — 1 complex parameters tj’s in T(z) are uniquely determined by the condition that 
S(z ) — T(z) is a polynomial of degree to. Then we choose a € C and b £ C by hand to maximize 
the number of intersections between the two sets: 

r T = {z | ImT(z) = 0}, F s = {z | ReS(z) = 0}. 

These intersections are the roots of the equation p n (z) + q m {z) = 0 where 

Pn(z) = S(z) + T(z ), q m (z) = S(z) - T(z), 

since 

Pn{z) + q m (z) = 2ReS'(z) + 2 ilmT(z). 




Figure 3. Curves Tt (black) and Tg (red) for for n = 4,5 and m = 2. The 
shaded region is fl. For n = 4 (left) we choose a = 0 and b = 1.1—0.1* to produce 
12 roots. For n = 5 (right) we choose a = 1.5 — 0.5* and b = —0.05 + 0.92* to 
produce 15 roots. 

Theorem 3 can be stated equivalently as follows: 
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Theorem (Theorem 3). For a = 0 and for a generic choice of b G C, the equation p n + q m = 0 
defined in terms of S and T (as above) has at least m 2 + m + n roots. 

Proof. Choose a = 0 and let Tg be the set of rays emanating from the origin and extending 
to the infinity, i.e., {oo x e^ +k ^ m ^ n : k = 0,1, ■ • • , 2n — 1}. Note that Ty has 2m + 2 curved 
rays emanating from b where every ray eventually approaches the infinity in the directions of 
{oo x e lm / n }i such that: i) different rays correspond to different values of l, and ii) l is chosen 
in a subset, that we will denote by N 2m + 2 ^ containing 2 m + 2 numbers from {0,1, • • • , 2n — 1}. 
Assuming that b Tg and that none of the rays hit any critical point of T(z) except at z = b, 
those curved rays do not intersect each other. 

Let Wf s (l = 0,1, • ■ • , 2n — 1) denote the connected components (that we will call “sectors”) 
in C \ Tg that contain the asymptotic direction oo x e Z7r */ n . A simple geometric consideration 
tells us that the number of intersections between Tg and a curved ray starting from b G W'V, and 
continuing into W{ 2 without passing the origin is at least 

min(|Z 2 — h\,2n~ \l 2 - /i|). 

This means that the minimal possible number of intersections between Tg and “the 2m + 2 
curved rays in T-t” is 

Y, min (l ^2 — Zi|,2n — \l 2 — Zi|) > 0 + 2(1 + 2 4-f m) + (m + 1) = (m + l) 2 . 

h€:N2m+2 

The remaining part of IV (i.e. that is not connected to b) approaches the infinity in 2 n — 2m — 2 
different sectors among Wi , that are not already taken by the rays from b. Assuming that there 
are no critical points of T in IV except the one at b , each curve in the “remaining part of IV” 
must connect two sectors from the 2 n— 2m — 2 sectors such that, around oo, each sector is ’’hit” 
by one curve only. Since there is at least one intersection between each curve and Tg, the minimal 
number of intersections between the “remaining part of IV” and Tg is n — m— 1 and the minimal 
number of intersections between IV and Fg is given by (m + l) 2 + n — m — 1 = to 2 + m + n. ■ 



Figure 4. Roots and II (shaded region) for (n,m) = (9, 7), a = 0 and b = 1.1 — 0.1L 
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2.4. A remark on Wilmshurst’s conjecture. Comparing with Wilmshurst’s conjecture, The¬ 
orem 3 undercounts the number of roots by 

(2) 3n — 2 + m(m — 1) — ( m 2 +m + n) = 2 (n — m — 1). 


One can show that the corresponding lemniscate has 2 m + 2 curves that connects b and a 
(cf. Figure 4). The numerics suggests that the 2 in + 1 regions in between these curves have, 
respectively, 

m — 1 , m — 2 , • • • , 2 , 1 , 1 , 2 , • • • , m — 1 , m, 


(counting from top to bottom in Figure 4) roots. Among these exactly 1+2+- • - Am = m(m+l)/2 
of them are found in the sense-reversing region (with to components; the shaded part in Figure 
4) and, therefore, the total number of roots must be at least 


( 3 ) 


2 x 


TO' 


+ m 2 

-b n = to 


+ to + n 


by Fact 1, giving the same number as in Theorem 3. 

Since there can be at most n— 1 components in 12, there can be n — to — 1 extra components in 
the sense-reversing region, and these components are not connected to the point b. Wilmshurst’s 
count is obtained when each of these components has exactly one root. This increases the total 
number of sense-reversing roots by n — to — 1 and, hence, increases the total number of roots by 
2 (n — to — 1), cf. (3), which is precisely Wilmshurst’s count, cf. (2). 

The various counterexamples studied in [9, 6] indicate that, the n — to — 1 “extra components” 
of 12 (that are not connected to b) can have more than one root in each component. For example, 
Figure 4 shows that there are two roots inside the component of 12 that is not connected to b. 

For to = 7i — 2, our numerical experiment supports the structure shown in Figure 4: there are 
to 2 + to + 7i roots that are counted in terms of the 2 to + 2 curves connecting a and b, and the 
excessive zeros are twice the number of zeros found in the component of 12 that is not connected 
to b. Choosing a = 0 and b = e l ' K ^ 2n ' 1 + e (e ^ 0 is needed so that the origin is not a root), 
we found that the number of “excessive zeros”, i.e. (the total number of zeros) —(to 2 + m + n), 
increases by 4 whenever n hits the numbers: 


7,15,22,30,37,45,52,60,68,75,83,90,98,105,113,120,128,136, • 


For example, for n = 100, there are total 13 numbers before 100 from the list, and the total 
number of zeros is 4 x 13 + (to 2 + m + n) = 52 + 9998, exceeding Wilmshurst’s count by 
4 x 13 — 2(n — m — 1) = 52 — 2 = 50. These experiments prompt us to suggest the following 
conjecture. 


Conjecture. When m = n — 2, the maximal number of roots of h = p n + q m is given by 

2 3 / \ 
n — —n + o{n) 

as n grows to oo (which is larger than Wilmshurst’s count of n 2 — 2n + 4). 


3. Geometry of caustics: Proof of Theorem 4 
As before, let 12 be defined by (1). 

Lemma 4. Setting 

f(z) =Pn(z)/q' m (z), 

let D be a connected component of 12 with exactly k (counting the multiplicities) zeros of f in 
D. Then, f : D —> O, where D is the unit disk, is a branched covering of degree k. 
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Proof. By definition, D is a connected component of / _1 (B). By the argument principle, for 
any w £ D, the number of preimages of w under / inside D is given by the winding number, 
^Aqd arg(f(z) — w). Since f(dD) C <90 and cannot “backtrack” on dD, the winding number 
does not depend onu/SD and it is k at w = 0 because D contains exactly k zeros of /. ■ 


When D contains exactly one critical point of p n , f : D —> D is a univalent map. In this case, 
let r] : [0, 27r) —» dD be the parametrization of <9 D given by 

r 1 {6) = r 1 {e' e ), 

where f^ 1 on 90 is obtained by the continuous extension of /: D — > D. This parametrization 
of dD is given by the harmonic measure of D (normalized by the factor 27r) with respect to the 
pole at the zero of /, i.e., 

(4) dd = dargf(z), z £ dD. 


This viewpoint can be be generalized when there are k zeros of / in D. In this case, the 
same equation (4) defines the parametrization, rj : [0,27rfc) -A dD , by the harmonic measure of 
fc-sheeted disk. 

On a smooth part of the curve dD (i.e., where f ^ 0), let v(z) be the tangent vector of the 
curve dD at z £ dD given by 


( 5 ) 


v(z) 


dy{9) = . f(z) 
dd f(z)' 


Let us consider the image of dD under h. For each z £ dD with f'(z) ^ 0 we obtain a tangent 
vector of the curve h{dD) at h(z) as follows (v is the tangent vector to dD, cf. (5), |/| = 1 on 
dD ): 


(6) V(z) = (v{z)d + v{z) 9) h(z) = v(z)p' n (z) - v(z) q' m (z), 

assuming this expression does not vanish. The next two lemmas deal with the geometry of the 
caustic that we identify now. 


Definition. The image of the critical lemniscate h(dfl) is called the caustic. 

Lemma 5. Let D be a component of Q with exactly k zeros (counting multiplicity) of f and 
dD is parametrized hy 77 : [0, 27 t k) — > dD such that d arg / (?y( 0 )) = d6. At z = r](6) £ dD where 
f(z) yf 0 and lm(v(z)q' m (z)^/ f(z)) ^ 0, we have 

d arg V(t](6)) 1 

— de = 2 

In other words, the caustic, away from the possible singularities, has constant curvature with 
respect to the special parametrization defined above. 


Proof. Note that, for 2 £ dD , the two terms in the right hand side of (6) have the same modulus 
(i.e., \p' n \ = \q' m \ on dD). We may rewrite V as 


V(z) 






19rn W12* Im 



= 2i VTRlm (v{z)q' m {z)y/J(zfj . 


( 7 ) 
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(Note that the result is independent of the branch of the square root function as \J f(z) appears 
twice.) If lm(v(z)q' m (z)y/f(z)j ^ 0, then we have, modulo 7r, 

1 7T —% 7r 

argF(z) = -arg f(z) + - = — log f{z) + -, 

where, in the last equality, we used again that \f(z) \ = 1 at z € dD. We have: 

da,TgV(’ij(9)) —i f di){6 ) 1 

= YJ d9 = 2’ 


( 8 ) 


dd 

where we used (5) in the last equality. 


The next lemma characterizes the “possible singularities”. 

Lemma 6. The only singularities of the curve h(dD) are cusps. When Zq £ dD is not in the 
branch cut of \ff there is a cusp at h(zo) if and only if the mapping from dD to R. given by 

z lm(v{z)q' m {z)\/ f( z )) on dD , 

changes sign across zq. When zq is in the branch cut of yff, the cusp occurs at h(zo) if and only 
if the mapping, 

777) \ 

7/7 o) j 

changes sign across zq. 

Proof. First note that away from the branch cuts of yfj, the set 

dDn{z: lm(v(z)q' m (z)y/J(zj) = 0} 

is finite. Indeed, otherwise as a level set of a harmonic function it would contain an arc. By (6), 
V = 0 over that arc and, therefore, h maps the arc into a point. Yet h, a harmonic polynomial 
of degree n, has finite valence < n 2 , a contradiction. 

If /' 7 ^ 0 on dD then all the functions (i.e. f(z), v(z ) and q' m (z )) that appear in the expression 
(7) for the tangent vector, V , of the caustic are smooth along dD. Therefore, from (7), a 
singularity of h(dD ) occurs only when V changes sign, hence h{dD) must have a cusp there. 
The lemma follows immediately. 

If f'(zo) = 0 for zq £ dD (i.e., it is a critical point of the lemniscate), then |u(2o)| = oo by 
(5). However, the argument of V (i.e., argY) is still either continuous or jumps by i r. And a 
cusp occurs when V changes direction, i.e., arg V jumps by 7r. ■ 

If f'(zo) = 0 for Zq € dD, multiple components of H merge together at Zo- It turns out that 
the image of dfl under h has an interesting structure, as we will see below. 


lm(v(z)q' m (z)y/J{zfj on dD , 


x H> Re 


Remark (a critical point on the critical lemniscate). Let z 0 £ dfl satisfy f(zo) = 
/( fc - 1 )7o) = Oand/( fe )7o)^0, i.e., 

/7) = /7o) + 7 - ^o) fc + 0(7 - a7 +1 )- 

Taking the absolute value and rotating if necessary, we obtain 


(9) 


1/7)1 = l/7o)l 


l/7o)l 

k\ 


Re 


/ (fc) 7o) 

/7o) 


{z - Z 0 ) k 


+ O((z~z 0 ) k+1 ). 


Locally, the lemniscate consists of 2k curved rays meeting at Zq and it divides the plane into 2k 
wedge-shaped sections with the angle ir/k at Zq. Among them, total k sections, equally spaced, 
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Figure 5. Lemniscate (left) and caustics, h(z) 


4-z- e ie z for 0 = 0 and 0 = £. 

4 6 


are included in Cl. To be more precise, taking a sufficiently small disk B centered at z$, B n Cl 
has exactly k components such that B D dfl consists of 2k curves emanating from Zq with the 
angular directions given by 


9j = r arg 


/0 o) 


k fW(zo) k 


+ T I 3 + l 


3 = 0 , 1 , 


,2k — 1. 


Note that the sections between Oil and Q-il+i for every l = 0,1, • • • , k — 1 are in Cl. Therefore, 
taking the single component (let us denote it by Dl ) of B n Cl between 0 2 e and 6 2 l+i, the 
tangent vector of its boundary changes angular direction from 0 2 l+i — it to 0 2 l at zq. Using (7), 
assuming zq is not in the branch cut of yff, the corresponding tangent vector of h(dDi) changes 
the direction by 7r at h(zo) (i.e., has a cusp singularity) if and only if 

(e l(e2<+ i ~ 7r) <L(z 0 )\/f( z oj) Im {e l9 ™q' m (z 0 )^/ /(z 0 )) < 0 


Im I 


or, equivalently, 


Oil + arg (q' m (z 0 )y/ f(zoj) £ (V ) u ( 


2k- 1 

7T, -;-7T 


When k is even (respectively odd) there can be at most two (respectively one) values of Us that 
do not satisfy the above condition. It means that for those values of Us h{dDi) has a smooth 
boundary at h(z 0 ), and for the other values of Us, h(dD e ) has a cusp at h(z 0 ). In Figure 5, the 
middle picture shows a caustic where the critical point Zq = 0 corresponds to the three cusps 
at the origin, and the last picture shows a caustic where one component (red) of the lemniscate 
maps zq to a regular point of the caustic. 


Lemma 7. Let D be a simply connected component of Cl with exactly k zeros of f. The number 
of cusps in h(dD) is odd (resp. even) when k is odd (resp. even) and, moreover, is >2 + k. 

Proof. According to Lemma 6 a cusp in h(dD) occurs whenever Im (vq' m y/J) changes sign. To 
locate such events, it is convenient to define the function, U/ : [0, 2nk) —>• M by 

(10) T : f? —>• arg fv(z) q m (z) vTw) . , = arg v(q(0)) + arg q' m (v{9)) + arg y/f(i](9)). 

For the term arg q' m (r/(9)), we choose the branch of the function “arg” such that the term is 
continuous with respect to 0. For the term yjf(rj(9)) we choose the branch of v7 and “arg” 
function (separately from the previous one), so that that term is continuous. Lastly, for the 
first term, argu, we choose the branch of “arg” such that the term is a piecewise continuous 
function where the only discontinuities are at the critical points of the lemniscate (i.e., where 
f(v(0)) =0). At the discontinuity arg v jumps by a positive angle, £ [7r/2,7r) from 9 2 l+i — tt 
to 0 2 l using the notations in the above Remark “a critical point on the critical lemniscate”; the 
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angle can vary according to the order k > 2 of the critical point. As a consequence, <F is 
a continuous function with only discontinuities being the jump(s) by angles in [tt/ 2,7r] at the 
critical points of the lemniscate. 

For a cusp to occur at 9, the condition in Lemma 6 gives that 

(11) (\F(0 — e), \F(0 + e)) contains In for some integer l. 

e>0 

Let us look at the three terms in if individually. The last term is linear in 9 with the slope 
1/2 by (8) and, therefore, Aargy/f = kn. (Here and below A stands for the increment over 
dD.) The second term, argq^, is a continuous function and A&Tgq' m = 0 because q' m does not 
vanish in the closure of D. Lastly, the first term, argv, is a piecewise continuous function with 
A arg v = 2n where the only discontinuities are at the critical points of the lemniscate (i.e., where 
f(v(9)) = 0). Summing up, the total increment of *F over [0, 2fc7r) is (2 + k)n and, therefore, 
there are at least k + 2 points (and exactly k + 2 if T is monotone) where the condition (11) 
holds. 

Suppose now that ’F is not monotone. For any point 9 that satisfies the condition (11) there 
are two possibilities: ase-> +0, 

(A): \F(0 — e) < In and + e) > In or (B): 'S’(9 — e) > In and ^(9 + e) < In. 

Since the total increment of *F is (fc + 2)-7r, the number of points satisfying the condition (A) must 
be larger than those satisfying (B) by exactly k + 2. Therefore the number of points satisfying 
either (A) or (B) can be bigger than k + 2 by an even number. As a consequence, the number 
of cusps is always odd (respectively even) when k is odd (respectively even). ■ 

Lemma 8. Let D be a component of Q with a single zero of f. If h(dD) has only three cusps, 
then it is a Jordan curve. There are more than three cusps in h(dD) (hence five or more according 
to Lemma 1) if and only if there exists a point with the winding number of h(dD) bigger than 
one, i.e., there exists p £ C such that A^d arg(/i(-) — p) > 2n. 

Proof. To prove the first statement of the lemma, it suffices to show that the three smooth 
(open) arcs between the three cusps do not intersect each other. Choosing any two arcs, there 
exists a cusp where the two arcs meet. Since the tangent vector rotates at most by 7r over a 
smooth part of the curve (cf. Lemma 5), the two arcs can only get farther from each other as 
one moves along the arcs starting from the common cusp. (Note however that, if the constant 
curvature is bigger than 1/2, h(dD) can self-intersect.) 

To prove the second statement, assume h(dD) has five or more cusps. Let T be a smooth curve 
that is obtained by slightly “smoothing” all the cusps of h(dD ), see the left picture in Figure 
6. Considering infinitessimally small smoothing, such deformation indicates that the tangent 
vector of T rotates by 

- • 2-7T — /£{cusps}7r < n — 5n = —4-71 

over the whole curve. Therefore T cannot be a (smooth) Jordan curve and h(dD) must have a 
self-intersection (that cannot be removed by a small perturbation, such as the one in the right 
picture of Figure 6). 

Let us define the orientation on h(dD) by the orientation inherited from dD. For each point 
p h{dD), one can consider the winding number of h(dD) around p. The winding number of 
h(dD) on the left side (with respect to the orientation) is bigger than the one on the right side 
by +1. Then, in a neighborhood of a self-intersection, there must be a pair of regions where the 
winding numbers differ by 2, see the right picture in Figure 6. ■ 
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Figure 6 . Cusp in caustic (left). Dotted line shows a “cusp after smoothing”. 
Relative winding numbers in a region around an intersection point (right). 


Now we prove Theorem 4. 

To have more than three cusps, the function 4/ that is defined in the proof of Lemma 7 must 
be non-monotonic, i.e., there must be a point where the slope of the graph of 41 is negative, i.e. 

nro d^(9) _ daxgv(r}(6)) , daxgq' m (r){0)) , 1 

[ 1 dO ~ d6 + dO + 2 < U ' 

Note that the first two terms in the left hand side give the curvature of q m (dD) with respect 
to the parametrization by 9 or, equivalently, re/|/'| where k is the curvature with respect to the 
arclength on dD (cf. (5)): 


dv(O) = d_ d qrn (q{9)) = darg 11 ( 77 ( 6 *)) daigq' m (q(0)) 

dd K de arg de de do 


So the inequality (12) is exactly the one that appears in i) of Theorem 4. 

Once the graph of 4/ is non-monotonic, one can create more roots of = 0 mod 7 r by 
vertically shifting this graph. This is done by the following transformation, which preserves the 
lemniscate {z : \p' n (z) / q' m {z)\ = 1 }: 

Pn e l<p Pn, P G R- 


Under this transformation, 

arg(^?Lv / 7) a*g(vq' m '/f) + V>/2 

because of the term ypf. This means that, for some ip, h(dD) can have more than five cusps. 
By Lemma 8 , there exists a point, say p £ C, where A dD(h{-) — p) > 2n. This means that 

Pn{z) = e lv p n (z) -p 

satisfies A g^h > 2n where h = p n +q^ and, therefore, h has at least two roots inside D. This 
ends the proof of Theorem 4. 


4. Construction of non-convex lemniscate 

Here we explain how we found the example in Figure 1. 

Theorem 5. Let f be a rational function and H = {z : \.f(z)\ < 1}. Let z 0 € dfl, f'(zo) = 0 
and (log fy'(zo) ^ 0. In a neighborhood of zq, dfl is a union of two smooth arcs that intersect 
perpendicularly at zq. Moreover, the point zq is not an inflection point (i.e. the curvature is 
strictly positive or negative at Zq) for either of the arcs if and only if 

(+in/4 (log fY"{zo) \ 

V (iog/)"U) 3/2 y 


Re 


^ 0 . 
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Figure 7. Lemniscate near zq, see text below Theorem 5. 


Corollary 2. If Zq £ 1 satisfies the assumptions in Theorem 5 and, furthermore, is not an 

inflection point ofdfl, then there is a connected component offl whose curvature of the boundary 
(with respect to the counterclockwise orientation) converges to a negative value at Zq. 

Proof. According to Theorem 5 there exists an open neighborhood U of zq such that SI D U is 
the disjoint union of two domains, see Figure 7 for an illustration (fl D U is the shaded region). 
When the “no inflection” condition in Theorem 5 is satisfied, the local configuration of fl has 
two possibilities: QC\U is the disjoint union of a convex and a non-convex domain (left in Figure 
7) or, S7 n U is the disjoint union of two non-convex domains (right in Figure 7). In either case, 
there exists a connected component of Q D U whose boundary has a negative curvature in a 
neighborhood of z 0 . ■ 

Remark. To relate Theorem 5 to our harmonic polynomials, we take f(z) = p' n (z)/q' m (z). For 
m = 1, as z approaches z 0 along the “concave” arc, k(z) < 0 by Corollary 2 and \f'(z)\ —> 0. 
Therefore, one obtains that k{z) /\f (z)\ —>■ —oo as z approaches zo along the “concave” arc. 
Hence we obtain the “non-convex” domain D as is needed to satisfy the condition n/\f'\ < —1/2 
in Theorem 4. 


In the rest of this section, we prove Theorem 5. 

Let / be a meromorphic function. Let 7 : K. —► C be the parametrization of a curve in C such 
that 7 ( 0 ) = zq and f'(zo) = 0. Taking the derivatives of log|/( 7 (f))| at t = 0 we obtain: 


^l°g|/(7(*))l 


= 0 , 

t =0 


^log|/( 7 (t))| 


^log|/(7(t))l 


t =0 


= Re(7(0) 2 (log/)"(z 0 )), 

He (^ 37 ( 0 ) 7 ( 0 ) (log f)"(zo) + 7 ( 0 ) 3 (log f)'''(zo) 


Above, 7 and 7 denote, respectively, the derivative and the second derivative of 7 with respect 
to t. Consider the Taylor expansion of log |/( 7 (t))| around t = 0. We have 

f 2 

¥ 


log|/(7W)l=Re(7(0) 2 (log/)"(z 0 ))^ 

+ Re( 37 ( 0 ) 7 ( 0 )(log /)"(z 0 ) + 7 ( 0 ) 3 (log fY"(z 0 )) | + 0(f 4 ). 


If 7 parmametrizes the lemniscate dfl passing through zq, where Q is as in (1), then we have 
log|/(7(f))| = 0 identically for all t, and all the coefficients in the above Taylor series must 
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vanish. The first coefficient vanishes when 

(13) 7 ( 0) 2 = fci(log/)"(z 0 ), ci G R. 

We may assume that 7 is the arclength parametrization, i.e., I 7 I = 1. Then 77 is purely 
imaginary and 

(14) 7 ( 0 ) = *c 2 7 ( 0 ) 


for some real constant C 2 . The second coefficient in the Taylor series gets simplified as 


Re 


(37(0)7(0)(log f)"{z 0 ) + 7(0) 3 (log /)'"(£ 0 )) 
- 3 cic 2 |(log/)"( 7 :o )| 2 - cilm(7(0)(log/)"(z 0 )(log/)" , (2:o)), 


using (13) and (14). The above expression vanishing implies that 


c 2 = --Im( 7 ( 0 ) 


Summarizing, we obtain 


(iog/r(zo) 

(log/)" ( 2 : 0 ) 

(log /)'"(*<>) 


7 ( ( ) = ,o + 7(0)<-^Plm^(0) (log/) „ (zo) 


t 2 + 0(t 3 ), 


where, using (13), 


MM - i V(lQg f)"(zo) _ ^ iv/i . \V0ogf)"(Z0)\ 

V(log/)"U)| V(!og /)"(*>) ' 

The two signs can be chosen arbitrarily and independently. This proves that there are two 
different arcs orthogonal at Zq. The curve d£l does not have an inflection point at zq when the 
quadratic term in the Taylor expansion of 7 (t) is non-zero, i.e., 


The proof is now complete. 


Re 


^ e ±W 4 


(log//"(so) \ 
(log /) ,, (~o) 3/2 / 


7 ^ 11 - 
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